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Abstract 



Charm photoproduction rates off nuclei below the nucleon threshold are estimated using 
the phenomenologically known structure functions both for x < 1 and x > 1. The rates 
rapidly fall below the threshold from values ~1 nb for Pb at the threshold (7.8 GeV) to ~lpb 
at 6 GeV. 

1 Introduction 

In view of the envisaged upgrade of the CEBAF facility up to 12 GeV it becomes important 
to have relatively secure predictions about the production rates of charm on nuclear targets 
below and immediately above the threshold for the nucleon target. This note aims at such 
predictions. From the start it has to be stressed that the dynamical picture of charm pro- 
duction at energies close to the threshold is much more complicated than at high energies 
(see e.g. [1]). Correspondingly the production rates cannot be described by the standard 
collinear factorization expression but involve gluon distributions both in x and transverse 
momentum, which are unknown at small transverse momenta (in the confinement region). 
To obtain our predictions we shall recur to some crude assumptions about these distributions 
both at x < 1 and x > 1. In the latter case we shall exploit the (poorly) known behaviour of 
the nuclear structure functions, which fall exponentially with x, the slope being in the range 
15-^16. With all these simplifications, we hope to be able to predict the rates up to factor 
2H-3. 



2 Kinematics 

Consider the exclusive process shown in Fig. 1 

-f + A ^ CC + A* (1) 

where A is the target nucleus of mass niA and A* is the recoil nuclear system of mass m\. 
The notation for all the momenta involved can be read off the Fig. 1. We denote the total 
mass of the CC system as M. Obviously M > 2mc where rric is the mass of the C-quark, 
which for our modest purpose we take as 1.5 GeV. The inclusive cross-section for charm 
photoproduction is obtained after summing over all states of the recoil nuclear system. 
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We choose a reference system in which the target nucleus is at rest and the incoming 
photon is moving along the z-axis in the opposite direction, so that q+ = q± = The (part 
of) the photoproduction cross-section corresponding to (1) is then obtained via the imaginary- 
part of the diagram in Fig. 1 as 

aA = Af 4^6{iAp - kf - mi')x \ ^ ) ag{M\e) (2) 




(27r)2 

Here V{k'^) is the vertex for gluon emission from the target; ag{M^,k'^) is the photoproduction 
cross-section off the virtual gluon of momentum k. We have also introduced the scaling 
variable for the gluon as x = k^/p^. Due to g+ = this is also the scaling variable for the 
observed charm. Note that this definition, which is standard at large eneries and produced 
masses, is not at all standard at moderate scales. In particular this x does not go to unity 
at the threshold for the nucleon target. Rather the limits for its variation converge to a 
common value 0.76. For the nuclear targets with A >> 1 its minimal value at the nucleon 
threshold is well below unity (-^ 0.64). One should have this in mind when associating this x 
with the gluon distribution: it follows that for a nuclear target, for energies going noticeably 
below the nucleon threshold, the contributing values of x still remain below x = 1. This 
eff'ect is essential for the charm cross-sections to remain not too small immediately below the 
threshold. 

We use the (5-function to integrate over k- to obtain imaginary part of the diagram in 
Fig. 1 as 

dxcfk^ fT{k^^ " 
2(A-x)(27r)3 I 



-A=A ,,r':^,,3 -9iM\ k-) (3) 



In these variables we find 
and 



= xsi A;^, si = 2pq (4) 



k^ = xAm- - - -^fci (5) 

where we have put j? = 'm}^ the nucleon mass squared, neglecting the binding. 

The limits of integration in (3) are determined by the condition > 4m^ , which leads 

to 

x{sx + An?) - -^—m*A - k\ - Am\ > (6) 

Since k\>0 one gets 

x{si + An?) - ^ m*y^ - 4m^ > (7) 

-/i X 

from which one finds the limits of integration in x: 

Xl < X < X2 (8) 

where 

xi^2 = —{As- ni*A + 4.ml ±^J[As- (m* + 2nic)^][As - {m* - 2nic)'^] (9) 

where s = si + An?. The limits of integration in k± at a given x are determined by (6). 
Using (5) we may pass from the integration variable kj_ to 

rxs,-4ml (^|jt2| fr{k'^)Y 1,2, ,2> 



Xl 



aA = i^ I xdx I „ 2 CFg{xsi -\k \,k ) (10) 



\k^m,u 2(27r)3 I P 
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where ^ 

\k \min — ~, X TT. (H) 
A — X 

The threshold energy corresponds to xi = X2 or = m\+2mc with the minimal possible 
value of the recoiling system, which is just vriA — Am. In terms of the photon energy E we 
have si = 2mE and the threshold energy is found to be 

It steadily falls with A from the nucleon target threshold 

Ef" = 2mc (^1 + 7.8 GeV (13) 

down to the value 2mc for infinitely heavy nucleus. 



3 High-energy limit 

To interprete the quantities cnctering Eq. (10) it is instructive to study its high-energy limit, 
which corresponds to taking si » ml and both quantities much greater than the nucleon 
mass. Assuming that the effective values of the gluon virtuality are limited (and small) one 
then gets for the nucleon target {A = 1) 




= nj^^ xdxj^ ^ j a,{xs,) (14) 

Here we also neglect the off-mass-shellness of the cross-section off the gluon, considering 
<< Am^. The obtained formula is precisely the standard collinear factorization formula 
with the identification 

xg{x,M^) = TT Z^*"' ||1 I^^J ^ j^^ d\k'\xp{x, \k'\) (15) 

Of course it is understood that this is only a part of the gluon distribution in x coming from 
a particular recoil state. The total gluon distribution is obtained after summing over all 
recoiling states. The quantity p{x, \k'^\) obviously has a meaning of the double distribution of 
gluons in x and Eq.(15) allows to relate function F^/A;^ entering our formulas of Section 
2 to the double distribution p(a;, \k'^\) and rewrite our formula for the cross-section at finite 
energies as 

aA = xdx / d\k'^\pA{x.,\k'^\)(Jg{xsi — \k'^\,k'^) (16) 

mzn 

4 The gluon distribution 

To make our estimates, we assume a simple factorizable form for the double density p{x, \k'^\) 
and choose the dependence in accordance to the perturbation theory with an infrared 
cutoff in the infrared region: 
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Function a{x) can be obtained matching (17) with the observed xg{x, M^) at a particular 
point Mq. Since we are interested in the threshold region, we take Mq = 2mc to finally obtain 

Our calculations show that the results are practically independent of the infrared cutoff for 
energies abov 6 GeV. However for lower energies they begin do rapidly fall with the growth 
of A. 

For the proton at x < 1 the gluon distribution g(x,Am'^) can be taken from numerous 
existing fits to the experimental data. In our calculations we have used GRV95 LO [2]. For 
the nucleus at x < 1 we use the simplest assumption gA{x,Q^) = Agi{x,Q^) neglecting the 
EMC effect in the first approximation. Obviously this not a very satisfactory approximation 
at X quite close to unity. 

For the nuclei at x > 1 we the gluon distribution may be estimated using, first, the 
existing data for the structure function in this region and, second, the popular hypothesis 
that at sufficiently low = sq the gluon distribution vanishes and the hadrons become 
constructed exclusively of quarks. Then one can, in principle, find the gluon distribution 
at a given from the standard DGLAP evolution equation with the quark distributions 
determined from the experimental data on the structure functions at x > 1 evolved back 
to = sq. In the present estimates we recurred to a simpler approach, neglecting the 
dependence of the quark distributions at a; > 1 in accordance with the naive parton model 
and experimental evidence for the structure functions at x > 1 which depend on very 
weakly. From [3] one concludes that at x > 1 the structure function of carbon falls with x 
as approximately exp(— Ax) with A ~ 16, its value at x = 1.05 being r-^ 6 * 10'^ and very 
weakly dependent on Q^. From these data one can find the quark density at x > 1 and using 
the DGLAP equation find the corresponding gluon density at x > 1: 

xgA{x,Q') = a^A'+'-'^-e-^-lnQho (19) 
ottA X 

with 5 = 16., a = 100. and sq the scale at which the gluon contents of the hadrons van- 
ishes. We take sq =(0.4 GeV)^. The ^-dependence has been chosen in accordance with the 
experimental data for hadron production at x > 1 [4]. 

5 Numerical results 

To simplify calculations we take in (16) the photon-gluon fusion cross-section on the gluon 
mass shell where it is known to be [5] 

where ec is the quark charge in units e, « = — — t and the limits t\^2 are given by 

ti,2 = -y[M±^M2-4m2] (21) 

We take the strong coupling constant ag = 0.3, infrared cutoff in the gluon distribution 
A = 0.3 GeV. 

With these parameters we obtain the cross-sections for charm photoproduction shown in 
Fig. 2. To have the idea of the number of nucleons which have to interact together to produce 
charm at fixed energy below threshold we show the limits of intergation xi and X2 in Fig. 3. 
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As expected the cross-sections rapidly fall for energies below threshold. Their energy 
dependence cannot be fit with a simple exponential (in fact they fall faster than the expo- 
nential). As to the aboslute values, for Pb at wc have obtained values around 1 nb at the 
nucleon threshold (7.8 GeV) and around 1 pb at 6 GeV. The A dependence is close to linear. 

6 Discussion 

We have crudely estimated charm photoproduction rates for nuclear targets below and in 
the vicinity of the nucleon target threshold. The estimates require knowledge of the gluon 
distribution in both x and k'j_ in a wide region of the momenta including the confinement 
region. 

Our estimates were based on a simple factorization assumption and introduction of an 
infrared cutoff. As mentioned at energies above 6 GeV the obtained cross-sections practically 
do not change at all when the cutoff is raised from 0.3 GeV to 0.4 GeV. However at 5 GeV 
they fall by nearly 30 times. Thus at energies much below the nucleon threshold and closer 
to the nuclear thresholds the results become strongly dependent on the cutoff. However at 
such energies the cross-sections becomes so small that do not present any interest from the 
practical point of view. 

The behaviour in the nearest vicinity of the threshold depends crucially on the choice of 
the scaling variable x and the form of the nuclear gluon distribution at x smaller but very 
close to unity. Our choice of x as the part of the "-|-" component of the hadron momentum 
carried by the gluon is well founded from the theoretical point of view. However with this 
definition x remains below unity at the nucleon target threshold. This considerably enhances 
the cross-sections for nuclear targets in the vicinity and somewhat below the threshold. Our 
assumption that in this region the nuclear structure function is roughly proportional to A 
is quite crude without any doubt. To improve the estimates one has to specially investigate 
this region of x for a nuclear target following [1]. 
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9 Figure captions 

Fig. 1. The forward scattering amplitude corresponding to reaction (1). 

Fig. 2. The charm photoproduction cross-sections for different photon energies E and 
targets. Curves from bottom to top correspond to A =1, 12, 64 and 207. 

Fig. 3. The limits of x-integration for different photon energies and nuclear targets. 
Curves from bottom to top correspond to A =12, 64 and 207. 




Fig. 1 



100 r 



CO 




EGeV 



1 / / 


1 1 1 1 1 1 


1 


1 




/ / 










' / 










,' / 










/' / 










/' / 










/ / 










/ 








- 


/ / 










/ / 










/ / 










/ 










/ 










/ / 










/ / 
_ ,' / 










'' / 








- 


; / 
; / 
,' / 










' ' 

; ■ 
1 . 

\ \ 








— 


\ \ i 
\\ '■ 










\\ 










\\ 










K\ 










— \ 








- 


" •X 

" - — _ * - 




















1 


1 1 1 1 1 1 


1 


1 





3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 

EGeV 



